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Abstract: We calculate the CMB temperature distortion due to the energetic electrons 



^ I and positrons produced by dark matter (DM) annihilation (Sunyaev-Zel'dovich effect, 

-vj ' SZdm) in dwarf spheroidal galaxies (dSphs). In the calculation we have included two 

cn ■ important effects which were previously ignored. First we show that the e with energy 

. ■ less than ~ GeV, which were neglected in previous calculation, could contribute a signif- 

^— V . icant fraction of the total signal. Secondly we also consider the full effects of diffusion 

^D . loss, which could significantly reduce the density of e at the center of cuspy halos. For 

neutralinos, we find that detecting such kind of SZ effect is beyond the capability of the 
current or even the next generation experiments, which is consistent with the quantitative 
/\ ' description made by S. Colafrancesco in [y. In the case of light dark matter (LDM) the 

C^ ■ signal is much larger, but even in this case it is only marginally detectable with the next 

generation of experiment such as ALMA. We conclude that similar to the case of galaxy 
clusters, in the dwarf galaxies the SZdm effect is not a strong probe of DM annihilations. 
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1. Introduction 

The physical nature of dark matter (DM) is a great unsolved problem in modern cosmology. 
A major observational approach to this problem is to look for possible signatures of DM 
annihilations. In many theoretical models, the DM particles could annihilate and produce 
7-ray photons and other energetic particles such as electrons and positrons. This is the 
case, for example, for the supersymmetric DM model [Q], and the light dark matter (LDM) 
model [^. For some early investigations on this subject, see e.g. [^, ^, ^, ^, ^]. 

Here we consider the energetic electrons and positrons produced by such annihila- 
tions. One possible way to reveal or constrain the presence of such energetic electrons 
and positrons is to look for the Sunyaev-Zel'dovich (SZ) effect induced by them (SZdm)- 
Energetic electrons or positrons could scatter with cosmic microwave background (CMB) 
photons, producing a distortion in the CMB spectrum. In clusters of galaxies, the virial 
temperature is sufficiently high that this effect could be produced by thermal electrons, 
which is the so called thermal SZ effect j^. The bulk movement of electrons with respect to 
the CMB also cause a similar effect, which is dubbed the kinetic SZ effect [|^]. In galaxies, 
however, the virial temperature is too low for the thermal electrons to produce significant 
thermal SZ effect, and the number of free electrons is also too low to produce significant 
kinetic SZ effect. However, non-thermal energetic electrons and positrons do exist such 
as cosmic rays in galaxies, which may produce an additional non-thermal SZ effect. In 
particular, DM annihilation may supply some of such energetic particles. However, DM 
annihilation is not necessarily the only source of cosmic ray particles. A number of "nor- 
mal" astrophysical sites, such as supernova remnants, AGN jets, pulsar winds etc., may 



also be partially responsible for cosmic rays. This is why the dwarf spheroidal galaxies 
(dSphs) may be of particular interest to us, as they are DM dominated and also considered 
to be very inactive, hence observations are less affected by background originating from 
normal astrophysical processes such as star formation and AGN activity. 

The SZ effect induced by a non-thermal distribution of charged particles were studied 
in 10, 12|, and then applied to clusters of galaxies [13, 14, 15] and dSphs, where the Draco 



dwarf was taken as an prime example [16|. These studies suggest that the DM- induced SZ 



effect could be used to search for signatures of DM annihilation. However, in a companion 
paper [^ (hereafter Paper I), we calculated the DM-induced SZ effect for clusters. We 
showed that contrary to previous claims, the DM annihilation induced SZ effect in galaxy 



clusters is not as large as previously work [16|, and there is little hope to detect it in the 



foreseeable future, which confirm the quantitative conclusion in ||]. Similar conclusion is 
also derived in a recent study [0]. In the present paper, we generalize our earlier work to 
the case of dSphs, for two DM candidates: the neutralino and the LDM. 

The rest of the paper is organized as follows: in Section 2, we discuss our method of 
calculation. We first briefly review the equations for SZ calculation with a given electron 
distribution. Then we discuss the density distribution of the DM in dSphs. Lastly the 
production of electron s and positrons from DM annihilation and propagation of these 
e^ are discussed. This differs somewhat from the case of clusters, as diffusion is more 
significant. In Section 3, we present the results of our numerical calculation. Finally we 
discuss our results and conclude in Section 4. 

2. Method of Calculation 

2.1 The nonthermal SZ Eflfect 

For the calculation of the SZ effect, our method here is the same as in Paper I (see also 



1 11]). The spectral distortion of CMB after travelling through a population of electrons or 

positrons ^ is given by 

Ar(x,g) {e- - If 

7^ = T1^9{x)y{0), 2.1 

Jo x^e^ 

where x = hu/kTo is the dimensionless frequency of CMB photon. To = 2.725 K is the 
undistorted CMB temperature, g{x) is the spectral distortion function, and y{9) is the 
Comptonization parameter for angle separation from the center. The spectral distortion 
function is given by 

9{x) = ^ [y ^o{xe-nPi{s)ds - ^o{x)^ , (2.2) 

where io{x) = x^/{e^ — 1) is the Planckian distribution of CMB photons, s = ln(z/'/i/) is 
the frequency shift of one photon after one scattering with e^. Here -Pi(s) is the frequency 
shift probability distribution after one scattering. 



Pl{s)= f fe{k)Ps{s,k)dk, 



^As all the calculations are the same for electrons and positrons, in this paper, the subscript e represents 
both e^ and e~ 



where k is the dimensionless momentum (momentum in units of mec), Ps{s,k) is the 
probabihty of a photon to have a frequency shift s when colhding with an electron with 
momentum k, fe{k) is the normahzed momentum spectrum of electrons 

1 dup dEp 



fe{k) 



rip dEp dk 



with ^^{Ep,r) the number density of e^ per unit energy interval, which will be discussed 
in details in the next section when the dSphs and DM candidate are specified. The number 
density of non-thermal electrons is 

^{E,,r)dEp, (2.3) 

with -Emax and -Emin the maximum energy and minimum energy of electrons produced by 
DM annihilation. The effective pressure of the non-thermal electrons can be defined as (see 
e.g. [H, |§) 



/■Kmax / k"^ \ 

^' = '"1 ''^'i^TTW )""''''• <'■'' 

with A^max = \/(-E'max/'"T-e)^ " 1> and fcmin = \/(-E'min/'"T'e)^ " 1- In the relativistic limit, 
this definition approaches to the usual approximation 

1 pl^m ax / \ 

Pe ^^{E), {E) = / fp{k) [Vl + k^ - l) mp dp, (2.5) 

in which (E) is the averaged kinetic energy. When DM annihilation is the only source of e^ , 
the upper limit of the integrals are set to be Smax = "i-dm and /cmax = A/("^DM/"^e)^ — 1- 



In some analyses fi^, Smin is set to be 0.01 ttt-dm, which was considered to be a good 
approximation as it is small enough comparing with the peak of the energy spectrum of 
continuous electrons (.Epeak = 1/20 moivi) as shown in ||^]. However, for the non-thermal SZ 
effect in dwarf galaxies, the contribution of electrons with E < O.OIttt-dm to the final signal 
is still significant. Therefore, we use Eq. (^^) which gives a more accurate expression for 
the contribution of trans-relativistic electrons and set -Emin = "T-e as no artificial cut off is 
necessary. Note however that even though this expression itself is accurate, there could still 
be significant error in the expression of feik) used for calculation at the lowest energies, 
see the discussion in §2.4 on e~^e~ production from DM annihilations. 

The contribution to Comptonization parameter by DM annihilation is given by the 
line-of-sight integral of effective gas pressure through the e cloud: 

y = ^[Pedl, (2.6) 

where ax is the Thomson cross section. 

The observed effective temperature change due to SZdm is an average of the temper- 
ature change within a beam 

ATobs{n,e) = J Bn{n',e)AT{n')dn', (2.7) 



where -Bn(n', 9) is the beam profile in the direction n' for a beam centered in the direction 
n with beam size 6. For illustration, we have made our calculation with a tophat beam 
profile and different beam widths. 

2.2 dark matter distribution 

In dSphs, DM dominates these systems both in the inner parts and outskirts. Similar as 
in Paper I, we consider the following three types of DM density profiles: 

p{r) = -^ — , , ,^, (hereafter B95, Ref.lp]), (2.8) 



P^'') = ( I ^n 1 / ^2 (hereafter NFW, Ref.@), (2.9) 



P^'') = ( I M.5hlr / U.51 (hereafter M99, Ref.[||). (2.10) 



(r/r,)i-5[l + (r/r,)i-5] 

They show similar behaviors (~ r^^) at large radii, but differs significantly near the center 
of the halo. We employ the B95, NFW and M99 profiles to represent the non-cuspy, 
moderately cuspy and strongly cuspy profiles of DM halos respectively. 

The density profile is truncated at a small radius where we assume the annihilating 



rate matches the in- falling rate of DM |24]. Within this radius, the density is kept at a 
constant value as 

IS "^nM 10~^^cm2s~^ q 

The density distribution of DM could in principle be determined from the velocity 
dispersion of stars. However, the uncertainty is large, and in the literature the required 
data is not always complete or easy to use. Parameters such as the characteristic density ps 
and radius r^ are chosen to be set by adopting specific density profile. Another generally 
used parameters are the virial mass Mvir and concentration parameters Cvir- The virial 
radius of a DM halo is defined as 

r -( Mv^^ V^' (2 12) 

where A is the overdensity, and pc ~ 139 Mq kpc~^ is the critical density of the universe. 
For the ACDM universe, A ^ ISvr^ + 82x - SQx"^ with 

X = 0,m{z) — 1 



is found to be a good approximation [^]. The concentration parameter Cvir is defined as 

Cvir = — , (2.13) 

where r_2 refers to the radius at which ^^^ ^ |r=r-2 = 0- The concentration parameter Cvir 
relates rvir and the density profile parameter as [gg] 

^B95 ^ ^vir ^NFW ^ ^ ^M99 ^ ^vir ^2 14) 



Simulations show that Cvir and Mvir are often correlated. Here we use a power-law relation 
between Cvir and M^ir as suggested in Ref. [27 

where cq = 9.6 and a = —0.1 for a ACDM cosmology. With this assumption, the DM 
distribution is fixed once Mvir is given. 

The halo of the dSphs would, however, be disrupted by the tidal force of the host 
galaxy (Milky Way), as a result part of the DM at the outer edge of the halo would be 
removed. To account for this effect, we calculate the tidal radius, which is given by the 
Roche criterion p8| 

MdSphjrt) MMw(D-rt) .„ ^ „. 

3 = — m V? — ' v2-16) 

rf (D - rt)'^ 

where D is the distance between the Galactic center and the dSph, rj is the tidal radius of 
the dSph and M{r) is the mass inside radius r. The results depend on the choice of the 



profile for Milky Way halo. Here we adopted the isothermal power law model [29] for the 
Milky Way 

^"° " 47rG (r2 + r2)2 ^^■^'> 

with Tg = 10 kpc and Va = 220 km s^^ as in |16, pO| . 

In Table ||, we list the parameters for a sample of 7 nearby dSphs, which were obtained 
by fitting the line-of-sight velocity dispersion profile of the stars [^]. NFW profile was 
assumed when deriving these parameters. Strictly speaking, these parameters should not 
be used for other halo profiles. To be self-consistent, one should re-estimate the virial 
mass with different halo profiles. However, since there are many uncertainties in both the 
observation and theoretical modeling, and the results of this paper should be regarded as 
an order-of-magnitude estimate, we will just use the same virial mass for M99 and B95 
profiles but re-evaluate the halo parameters. 

The nearby and relatively large Draco dwarf represents one of the best candidates for 
searching for DM in dSphs, and it has been studied extensively with many observations 
at different wavelengths. Since the basic properties of the SZ effects in different dSphs are 
qualitatively similar to each other, we will take the Draco dwarf as the prime example in 
the following discussions. 

2.3 e^ production from DM annihilation 

As the SZ effect produced by positrons and electrons can not be distinguished, here we use 
e to refer both e^ and e~ . The e source function from DM annihilation can be written as 

^ "'■DM f '^^e. 

where {av) is the velocity-weighted annihilation cross section, rriDM is the DM particle 
mass, p{r) is energy density of DM, and dN^/dE^ is the number of electrons produced per 



Table 1: The dSphs parameters used in this paper. NFW profile of dSph is assumed when 
calculating Vg, ps and rj. The distance and virial mass data are mostly taken from Ref. |3l[|, other 



parameters are calculated as outlined below. For Ursa Minor, which was not included in Ref. |31 
we take the mass as to be the same as that of Draco. 



Name 


D 

|kpc| 


Mvir 

flO»MQl 


ps 


[kpcl 


rt 


flO»M0/kpc^l 


[kpcl 


Draco 


80 


40 


0.82 


1.2 


9.9 


Leol 


250 


10 


1.2 


0.64 


16.7 


Fornax 


138 


10 


1.2 


0.64 


10.3 


LeoII 


205 


4 


1.5 


0.43 


10.6 


Carina 


101 


2 


1.8 


0.32 


4.8 


Sculptor 


79 


10 


1.2 


0.64 


6.5 


Sextans 


86 


3 


1.6 


0.38 


4.8 


Ursa Minor 


66 


40 


0.82 


1.2 


8.6 



annihilation per energy interval. In the following, we discuss two types of DM particles: 
neutralino in the supersymmetric model, and the LDM. 

Neutralino: We consider a neutralino with the typical values of parameters {av) = 
3 X 10^^ cm^ s~^ and m-DM = 100 GeV. The direct channel to e^e~ is generally suppressed 
for neutralino, so electrons are in most cases produced from the cascades of the annihilation 
final-state particles such as heavy leptons, quarks and gauge bosons j^. The spectra 
of electrons can be different from each other for different annihilation modes. We use 
the package DarkSUSY |32| to calculate the final-state spectra of electrons. Our fiducial 



annihilation mode is assumed to W~^W . 

In Fig. H, we plot the resulting electron energy spectrum for several annihilation modes. 



In the same figure, we also plot the parametrized fit of spectra used in |16| for comparison. 
As we can see from the figure, the parametrization fits the DarkSUSY spectrum well 
above 1 GeV. However below 0.1 GeV this analytical form deviates significantly from the 
DarkSUSY results. 

LDM: Inspired by the excess of 511 keV lines at the Galactic center observed by 



SPI/INTEGRAL H, g], it was proposed that a light DM particle with mass 1 - 100 MeV 



may be responsible for the data |35|. To allow the production of positrons through its 
annihilation, one needs m-DM > w-e- The mass of LDM should also be less than ~ lOOMeV 
to avoid producing too many 7-ray photons from n final state. The model has since been 



further constrained by a number of more detailed considerations on its phenomenology [36 



37| , P8| , |39| , [40| . For example, a small fraction of the energetic positrons produced in the 
annihilation could be directly annihilated into energetic 7-rays (in-flight annihilation), and 
the limit on 7-rays yields stringent limits of moM < 3MeV ||3^ or m-DM < 7.5MeV p9| , 
depending on the ionization state of the interstellar medium in the Galactic center. A 
comparable bound of m-DM < 7MeV was derived using CMB data by considering the effect 
of LDM annihilation during the recombination history [HO] . 

Here, as in Paper I, we assume the mass of the LDM is 5 MeV, which is consistent 




100 



Figure 1: Electron yield spectra dN/dE for W^W~ (solid curve), bb (dotted curve) and t^t^ (dot- 
dashed curve) annihilation modes for neutralino with mass m^ = 100 GeV. The dashed curve refers 
to the analytical fit adopted in |l6|. 



with the limits obtained in ||3^, |4^. It is slightly beyond the constraint of [^], but note 
that the SZ effect here would be stronger for smaller DM mass, so our choice is actually 
conservative. The cross section derived from the flux of the 511 keV line emission from 
the bright Galactic bulge region is (o-t;)LDM ~ 10^°(mDM/MeV)^ cm^ s^^. Thus for our 
mass choice of LDM we have (av) = 2.5 x 10~^^ cm^ s~^. The recent observations of the 



spatial morphology of the line emission by INTEGRAL |41] found a significant asymmetry 
in the disk emission with a resemblance to the observed distribution of low-mass X-ray 
binaries in the hard state, which indicates that these X-ray binaries might be the main 
sources. However, it is worthwhile to note that the large bulge-to-disk ratio (3 ~ 9) is 
not easily achieved even with the assumption that a large fraction of the disk positrons 
is transported via the regular magnetic field of the Galaxy into the bulge and annihilate 
there, so the LDM model can not be excluded yet. Our assumption about the annihilating 
cross section may be regarded as the upper limit of LDM contribution to the observed 511 
keV line. We assume the produced e^ have a monochrome spectra dNe/dEf, = (5(£' — ttidm) 
in this case. 

2.4 e propagation in dwarf galaxies 

DM annihilation injects electrons and positrons in the galaxy halo at a constant rate. The 
propagation of these charged particles in the tangled magnetic field can be modeled by 



diffusion. Tliey also lose energy by radiation during this process. As a result, the e 
spectrum satisfies the following transport equation [42|: 



d due 
didK 



V 



D{r,E)V 



dE, 



+ 



_d_ 

dE 



b{r, E) 



due 
dE, 



+ Qeir,E) , 



(2.19) 



where drie/dEg is the number density of e per unit energy interval, D{E, r) is the diffusion 
coefficient, b{E,r) = —dE^/dt represents the energy loss rate, and Qe{E,r) is the source 
function. For simplicity, we assume that D and b are independent of spatial location. The 
diffusion coefficient D{E) is assumed to have a power law dependence on energy E and 
magnetic field B: D{E) = Dq {E/B)\ with 5 = 1/3 |§, although it is not clear to what 
extent this relation could be extrapolated. Due to the much smaller scale of uniformity 
of the magnetic field in dwarf galaxies, we set Dq = 3.1 x 10^^ cm^ s~^, which is lower 
than the value in the Milky Way P, |4^. We also set B^ = 1 /iGs as our fiducial value of 
magnetic field. 

The energy loss rate is given by b{E) = biQ^[E) + bsyn{E) + bioniE), i.e. 



KE) 



10-i7GeV s 



- = '■''' (liSF)' + °-2=5^ (T^lil^)' + 2^i;^l'"(r) + 6.61, (2.20) 



where A'^h — 1-3 x 10 ^cm "^ is the number density of neutral gas in dSphs, T and /3 are 
the Lorentz factor and the velocity of e respectively. 
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Figure 2: Comparison of the time scales for various energy loss mechanisms and diffusion in a 
dSph galaxy with size r^ = 1.6 kpc. B = 1/iGs and Ah = 1.3 x lO^^cm^"^ are assumed. 



The characteristic time scales for different energy loss mechanisms and spatial diffusion 
are plotted in Fig. ^. The energy loss time-scale is defined as rioss = E/b{E), and the 



diffusion time-scale is r^if ~ rf^/D{E). We clearly see that tlie ionization loss has negligible 
effects, except below MeV energies. However, as the size of the dwarf galaxy is small, 
spatial diffusion could dominate at this energy range. Synchrotron emission and inverse 
Compton scattering (ICS) off the CMB photons are important for high energy electrons. 
The transport equation ( |2.19 ) could be solved with the Green's function method as 



described in 
form 



For time-independent source function, the equilibrium solution has the 
1 






m-DM 



dE' G{r,Av)Qeir,E'). 



(2.21) 



b{E) JE 

where Av{E,E') = J^ deD(e)/6(e) is the mean scale of diffusion covered by an electron 
while losing energy from E' to E. For a spherical symmetric system with free escape 
boundary condition at r/j (i.e. assuming the magnetic field is so weak that beyond this 
region the electrons are not confined), the Green's function is given by 



G{r,Av) = ^=^ y (-1)" / dr'— ^ 






exp 



(r 



.M2 



AAv 



exp 



(r -|- r: 



/^2 



4 At) 
(2.22) 



where r'^ = (— l)"r'-|-2nr/j is the location of nth "charge" image for r^ (c.f. Ref. |13]). The 
value of r^ is generally adopted as twice of the radius of the stellar component, typically a 
few kpc for local dSphs. 
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Figure 3: The Green's function G varies with respect to \/Av for two different radii r — 0.001 kpc 
and r = 0.1 kpc and for three density profiles: B95, NEW and M99. 

The behavior of the Green's function G (r. At;) is shown in Fig. y. We plot G as a 
function of \/Av for two different radii and for the three density profiles described in Sec. 



2.2. As expected, the Green's function decreases as the diffusion length VAv increases. 
Near the center of the galaxy, however, we find G <C 1 as represented by the r = 0.001 kpc 
curves for the M99 and NFW profiles, indicating that the diffusion effect is significant and 
it will affect the propagated spectrum. We can also see that for more cuspy density profile, 
the decrement in G is also larger. 

When spatial diffusion can be neglected, G ~ 1, we obtain the solution for no diffusion 
or in situ energy loss: 






nodif 



/ dE'Qe{r,E'), (2.23) 

Je 



HE) 



The spatial distribution of e then traces the source function, which is proportional to the 
density square of DM. However, for dSphs this approximation is generally invalid. 

3. Results 

3.1 Distribution of electrons 



To take the diffusion effect into account, Ref. |13] also proposed an analytic approximate 
solution of the transport equation 

^ (r, E) « [Qe{E, r)n,ss] x -^ x ^^ , (3.1) 

where V^ oc r| is the volume occupied by the DM source^, and Vo oc \^{E) is the volume 
occupied by a diffusing electron which travels a distance \{E) ss [D{E) x TiossiE)]^'"^ 
before losing a significant fraction of its initial energy. This changes the total number of 
electrons, but does not change the shape of its density profile. Here we call this model as 
the "approximate solution" . 

We plot in Fig. |^ the total number density distribution of e produced by neutralino 
annihilation for the case of dSph Draco with an assumed M99 profile. Our full solution 
obtained with the Green's function method is shown as the solid curve, while the "approx- 
imate solution" is shown as the dotted one. The dot-dashed curve below the solid one is 
the result obtained with i?min = O.Olm-DM- We see that the non-thermal electrons with 
E < O.OlmoM still make up a significant fraction of the total density. Previous analysis 
on SZdm in dwarf galaxies in Ref. |jl6| predicted /iK-level temperature distortion within 
1 arcsec for strongly cusped (M99) dark halo. They employed the approximation 



mDM ^n, 



"e(^) = / -—=;{r,E) 



O.OlmDM '^E 



AE. (3.2) 

nodif 



This "no-diffusion solution" is shown as the dashed line. It is apparent that there are sig- 
nificant differences between the full solution and the simplified solution or the approximate 
solution. The e density profile of our full solution has a significantly shallower slope than 



^As the galactic magnetic field has a limited extent, here we define ru as the radius of diffusion zone 
instead of the tidal radius of the dark halo. 
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Figure 4: The e^ density ne{r) from neutralino annihilations in the Draco dSph. The sohd and 
dot-dashed curves are our fuU solutions with -Emin = m-e and O.OIttidm respectively. The dotted 
curve is the result of the "approximate solution" , and the dashed curve is the result of the "no- 



diffusion solution" as described in Eq. (3.2) 



the "approximate solution" . Thus in the inner region, the full solution predicts much lower 
electron number density than either the "no-diffusion solution" or the "approximate solu- 
tion" . We can understand this by noting that diffusion could remove electrons in the inner 
region and spread them in a larger region. Although the "no-diffusion solution" attempts 
to account for the effect of diffusion by reducing the total number density of electrons, it 
reduces the electron density by the same factor at all radii. Actually the inner region will 
be affected much more significant than the outer regions. Also, for more cuspy halos, the 
scale over which p'^{r) has a significant gradient is smaller, resulting in more remarkable 
loss of electrons due to diffusion. 

In Fig. we plot the density distribution of electrons produced by DM annihilation 
for different DM halo density profiles, all of which are obtained with the Green's function 
method. As expected, the strongly cuspy M99 profile produces the largest e^ density, 
followed by NFW profile and then B95 profile. In the center, the density of electrons of 
M99 profile is several orders of magnitude higher than that of NFW profile. Thus, the 
density profile of DM halo is crucial in determining the strength of the SZdm signal. 

3.2 The SZdm signal 

The expected spectral distortion of CMB due to the energetic e^ from DM annihilation 
observed with 1" beam aimed at the center of the Draco dSph are shown in Fig. g. The left 
panel is for a 100 GeV neutralino, and right panel is for a 5 MeV LDM. The three curves 
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Figure 5: The e* density for different halo profiles assumed The red curve is the result for M99 
profile, blue for NFW and black for B95. 
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Figure 6: The SZ effect induced by the neutralino (left) and LDM (right) annihilations. The result 
for M99 profile (red), NFW profile (blue) and B95 profile (black) are shown. 

in the figure are for M99, NFW and B95 DM halo profiles respectively. As expected, 
the strongly cuspy profile (M99) yields greater temperature distortion, while NFW and 
B95 profiles yield curves which are similar but with less magnitude. We plotted only the 
absolute value of the temperature distortion in this figure, note that at low frequency (small 
x) we have AT < 0, while at high frequency (large x) AT > 0. As a result the magnitude 
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Figure 7: Spectral Distortion as a function of x = hiy/kTcMB for the case of Draco, with beam 
size ranging from 1" (red), 6" (green), 1' (blue), 6' (yellow) to 15' (black). 



of the temperature distortion appears to drop in the middle of the figure, where it reverses 
its sign. Unlike the case of the thermal SZ effect, this null point is not always fixed at 217 
GHz, but varying slightly. 

The effects of different beam sizes are illustrated in Fig. |^. The red, green, blue, yellow 
and black curves are the observed temperature distortion for the beam size 1", 6", 1', 6' 
and 15' respectively. The smaller beam probes a smaller region around the center, hence 
the SZdm signal is larger. 

We find that the magnitude of the signal is much smaller than the previous estimate in 
the literature (e.g. Ref.|jl6|). In Fig. ^ we plot the SZdm signal as a function of the beam 
width at X = 0.616 {v = 35GHz) which is similar to the corresponding value in Ref. ||16|. 
We also take the same values as TTiDM = lOGeV, < av >= 1.0 • 10~^^cm^ s~^ in order 
to have a direct comparison. The solid line is for DarkSUSY calculated electron injection 
spectrum with Green's function solution of the propagation equation. Other lines use the 
analytical form of the electron spectra which were shown as the dashed line in Fig. |l|. The 
dot-dashed line in Fig. |8| is the result from our full solution with the Green's function, 
while the dotted line is for the "approximate solution" , and the dashed line is for the "no- 
diffusion solution" respectively. Due to the large difference in the low energy electrons, the 
resulting SZdm signals are quite different between the DarkSUSY and the analytic spectra. 
On the other hand, if one neglects diffusion or only treats it as a global reduction in the 
total number of electrons within the halo as were done in previous works, the signal could 
be much higher near the center of the halos. With the diffusion effect included, however, 
the brightness temperature only increases slightly near the center. Only when the beam 
width become large enough, which is comparable with the region occupied by the stellar 
component, the approximate model can be used. For the 10 GeV neutralino, even for the 
M99 profile the expected signal in our fiducial model (shown as the solid line in Fig. ^ 
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Figure 8: SZqm signal as a function of observing beam width 9. The solid line is our fiducial 
model with DarkSUSY electron spectra. Other lines use analytical spectra, with different approach 
to solve the transport equation: Green's function solution (dot-dashed), "approximate solution" 
(dotted), and "no-diffusion solution" (dashed) respectively. 



at low frequency is well below 10~^/iK, which is not observable with the current or next 
generation instruments. At the high frequency end, the signal is higher, but it is difficult 
to observe from most sites on the ground. Furthermore, although the dSphs do not have 
a thermal SZ background, such small signal might be swamped by other foreground or 
background noises. 

Finally, we illustrate in Fig. |9|the effect of different diffusion parameters in D{E) which 
are not yet well-constrained. We plot the spectral distortion for three groups of Dq and 5. 
The solid line corresponds to our benchmark value Dq = 3.1 x 10^^ cm^ s~^ and 6 = 1/3. 
The dot-dashed line is for the same 6 but Dq = 1.0 x lO^^cm^ s^^. The dashed line adopts 
the same Dq but 6 = 2/3. For this purpose we only choose the LDM case to plot. We see 
that the SZdm effect decreases with increasing Dq or decreasing 6. But when the diffusion 
parameters vary within a reasonable range, the basic conclusions remain unchanged. 

In the above we have discussed various kinds of effects on SZdm, taking Draco for 
example. For other luminous dSphs as listed in Table |l], we present the values of AT at 
35 GHz and 1 THz in Table for the neutralino case, and in Table |3| for the LDM case. 
M99 profile is assumed. The particle physics parameters of DM are the same as described 
in Sec. 2.3. Since we have essentially assumed the same density profile for these dSphs, 
the magnitudes of their SZdm effects are determined mostly by the the mass Mvir of each 
galaxy. The differences in distance, on the other hand, do not affect the SZdm effect. Thus, 
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Figure 9: The temperature variation for beam size 1" induced by LDM annihilating in dark halo 
with M99 density profile. The solid line corresponds to our benchmark value Do = 3.1 x lO^^cm^ s~^ 
and (5 = 1/3. The dot-dashed line is for the same 8 but Dq — 1.0 x lO^^cm^ s^^, and the dashed 
line is for the same Dq but 6 = 2/3. . 



with these assumptions, dSphs with larger Myij. would generally induce larger temperature 
distortion. This can be easily seen by comparing the result for Draco and Leo I, which 
have the same r^ but different Mvir- Another essential factor in determining the final 
temperature distortion is the size of the diffusion zone in each dwarf galaxy. The value of 
r/j is generally adopted as twice of the radius of the stellar component (if this value is larger 
than the tidal radius, we set r^ = rt). For smaller rh, a larger fraction of e would escape 
away, and the SZqm effect is also smaller. This is reflected clearly through the comparison 
between Fornax and Leo L 



Recently, the PAMELA gj, ATIC g§ and FERMI |§ observations of the cosmic ray 
electron and positron energy spectrum indicate a possible excess or even peak in the cosmic 
ray electrons and positrons. If this excess is due to DM annihilations, then the annihilation 
cross section must be "boosted" by a large factor compared with the usual expectation. 
Such a boost could be due to, e.g., "Sommerfeld enhancement" |^, H, 49, ^, 51]. However, 



in such a case the mass of DM particle must be of TeV scale, so even with a 1000-fold 
enhancement, the induced SZdm effect, according to our estimate, would still be quite 
small. The DM substructures inside the dark halo are also thought to be able to "boost" 
the annihilation signal to some extent |52, |5^. However, the substructures near the center 
of the halo will be destroyed by tidal force. The "boost" is more effective at large radii 
1 54]. Furthermore, for the cuspy profiles the central density is high enough to dominate 
the contribution over that from substructures. Hence the "boost" , regardless of its origin. 
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Table 2: Neutralino-like DA^E induced SZ effect (in units of K) for Local Group luminous dSphs. 



We assume m^ = 100 GcV, (av) = 3.0 x 10" 


-26^j^3g-i^ M99 profile, 


and = 1". 


dSph 


rfe(kpc) 


M^UMq) 


35GHz {x = 0.616) 


lOOOGHz {x = 17.46) 


Ursa Minor 


1.6 


4x 10^ 


-5.44 X 10-12 


2.43 X 10-1° 


Draco" 


1.6 


4x 10^ 


-5.39 X 10-12 


2.32 X 10-10 


Leo I 


1.6 


1 X 10^ 


-6.88 X 10-13 


5.37 X 10-12 


Fornax 


5.4 


1 X 10^ 


-5.12 X 10-13 


1.91 X 10-1° 


Leo II 


1.04 


4x 10^ 


-5.45 X 10-15 


3.35 X 10-11 


Carina 


1.7 


2x 10^ 


-2.37 X 10-13 


8.97 X 10-13 


Sculptor 


3.26 


1 X 10^ 


-2.61 X 10-12 


9.72 X 10-11 


Sextans'* 


4.8 


3x 10^ 


-1.61 X 10-12 


5.89 X 10-11 



"^stellar IS around 0.93 kpc, which give a slightly large r/j as 1.86 kpc. In the calculation, 
we just use the r^ = 1.6 kpc as the diffusion zone for consistency with previous results, 
^rgteaier is around 4 kpc, which give a really large r^ as 8 kpc. In the real calculation, we 
use rt as the diffusion zone. 

Table 3: LDM induced SZ effect (in units of K) for Local Group luminous dSphs. We assume 
TOldm = 5 MeV, {av) = 2.5 x lO-^^cm^s"!, M99 profile, and 9 == 1". 



dSph 


r/,(kpc) 


M^UMq) 


35GHz {x = 0.616) 


lOOOGHz (x = 17.46) 


Ursa Minor 


1.6 


4x 10^ 


-3.11 X lO-'^ 


2.68 X 10-5 


Draco 


1.6 


4x 10^ 


-3.11 X lO-'^ 


2.67 X 10-5 


Leo I 


1.6 


1 X 10^ 


-6.76 X 10-^ 


1.48 X 10-^ 


Fornax 


5.4 


1 X 10^ 


-2.37 X lO-'^ 


3.80 X 10-5 


Leo II 


1.04 


4x IQS 


-1.63 X 10-^ 


4.76 X 10-6 


Carina 


1.7 


2 X 10^ 


-1.97 X 10-^ 


1.71 X 10-6 


Sculptor 


3.26 


1 X 10^ 


-1.37 X 10-^ 


1.48 X 10-5 


Sextans 


4.8 


3x 10^ 


-7.76 X 10-s 


1.09 X 10-5 



would have little impact on our present calculation. 

4. Conclusion 

In this paper, we calculate the non-thermal SZ effect induced by the energetic electrons 
and positrons produced by the annihilation of neutralino and LDM in dSphs. We take the 
Draco dwarf as an example to present most of our results, but we have also applied the 
same calculation to other luminous dSphs in the Local Group. 

In our calculation, we obtain the equilibrium e^ distribution by fully solving the dif- 
fusion equation with the Green's function method. We take synchrotron emission, ICS, 
ionization loss as well as diffusion loss of e into account. We find that for small scale 
systems such as the dwarf galaxies, diffusion effect is crucial. The e^ distribution is much 
less steep than the source function which is proportional to Pj^}^- This important effect 
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was mentioned in ||l| but was neglected in the numerical calculations in their previous anal- 
ysis 1 16 1, and the resulting e density could be different by several orders of magnitude. 
Another difference is that we considered contributions from e^ with relatively low energy 
(i.e. we do not use any arbitrary energy cut off). Such low energy e^ could still contribute 
significantly to the SZdm signal. 

The SZ effect by the non-thermal e depends on the density profile of the DM halo. We 
considered M99, NFW and B95 profiles. As expected, for the cuspy M99 profile the effect 
is much more significant. Nevertheless, for all profiles we find much smaller signals than 
previously claimed. Due to the diffusion effect, even for the strongly cuspy M99 profile, the 
SZ effect predicted for neutralino is too small to be observed with the current or coming 
generation of instruments. The possible astrophysical foregrounds or contamination may 
even make the detection more difficult. 

Considering that a smaller mass of DM will result in a higher DM number density and 
annihilation rate, we also investigated the case of LDM, i.e. DM candidate with MeV mass 
scale. In this case, there is some hope of detecting the DM induced SZ effect, though the 
effect is also much smaller than previous claims. For the strongly cusped profile assumed, 
AT could reach tens of ^K when the frequency is around or larger than THz. However, 
we note that THz observation is difficult to do on ground, with the possible exception of a 
few sites in the Antarctica high plateau. 

One of the largest uncertainty in our calculation is the diffusion coefficient of the 
dSphs, which depends on the magnetic field strength, for which we have no data. We made 
an assumption on its form and value for the present calculation. It is possible that the 
magnetic field is much weaker than we assumed, in that case, the SZ effect by non-thermal 
e would be even smaller. 

These results show that the non-thermal SZ effect induced by DM annihilation in dSphs 
is small and difficult to observe, even though the expected astrophysical contamination is 
small. Combined with our earlier results for clusters (Paper I), we conclude that the 
non-thermal SZdm effect is perhaps not a powerful method to detect or constrain DM 
annihilations. 
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